Gapless inhomogeneous superfluid phase with spin-dependent disorder 
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We show that the presence of a spin-dependent random potential in a superconductor or a super- 
fluid atomic gas leads to distinct transitions at which the energy gap and average order parameter 
vanish, generating an intermediate gapless superfluid phase, in marked contrast to the case of spin- 
symmetric randomness where no such gapless superfluid phase is seen. By allowing the pairing 
amplitude to become inhomogeneous, the gapless superconducting phase persists to considerably 
higher disorder compared with the prediction of Abrikosov-Gorkov. The low-lying excited states are 
located predominantly in regions where the pairing amplitude vanishes and coexist with the super- 
fluid regions with a finite pairing. Our results are based on inhomogeneous Bogoliubov-de Gennes 
mean field theory for a two dimensional attractive Hubbard model with spin-dependent disorder. 

PACS numbers: 74.62.En, 37.10.Jk, 74.40.Kb 



Introduction: The interplay of disorder and strong inter- 
actions and the mechanisms underlying metal-insulator 
transitions is still an open and very important question. 
Traditionally, this problem has been discussed in the con- 
text of solid state systems. More recently it has become 
possible to emulate quantum Hamiltonians using ultra- 
cold atomic gases [ij], opening up a new dialog between 
condensed matter and atomic physics leading to new in- 
sights into the effects of correlation on randomness. 

The Anderson model for localization [2] , first proposed 
in 1958, describes the possibility of localized electronic 
states formed by quantum interference of waves in a ran- 
dom potential. Later, in the famous "gang of 4" paper, it 
was shown using renormalization group flows of the con- 
ductance that a system in two or lower dimensions would 
always flow to a localized state for arbitrarily weak dis- 
order. In contrast, in three dimensions there is a critical 
disorder separating localized and conducting states 
Over the years the precise verification of the Anderson 
model has not been possible because in real systems in- 
teractions cannot be ignored, and, moreover, become in- 
creasingly important as the states get more localized near 
the metal-insulator transition. However, very recently 
in ultracold atomic gases the Anderson model has been 
emulated for the first time by tuning the interactions 
to zero [4] using Feshbach resonance and the localized 
wave function has been visualized. In these experiments, 
randomness is introduced through optical speckle [5] or 
bichromatic lattices [6]. 

If attractive interactions between fermions are turned 
on, such systems become superconducting (for charged 
electrons) or superfluid (for neutral fermions). These 
states can exist in two dimensions in a disordered system 
below a critical disorder even when all the single par- 
ticle states are localized. The ensuing superconductor- 
insulator transition has been studied in considerable de- 



tail in the literature j^-Q- A particularly intriguing as- 
pect of cold atoms is the ability to tune the interactions 
either through the Feshbach resonance or through opti- 
cal lattices. In the limit of strong attractive interactions 
between fermions the problem crosses over to bosons in 
a random potential lo|-l3 where new glassy phases are 
expected for large disorder. In this way, the complex 
interplay of interactions and disorder is combined with 
BGS to BEG crossover physics [15]. 

In this paper, we investigate the effect of a spin- 
dependent random potential on a superconductor. A 
related problem, that of magnetic impurities in a super- 
conductor, was studied by Abrikosov and Gorkov [l6| 
who found that the superfluid density was suppressed 
for fairly small impurity concentrations. These diagram- 
matic calculations were performed at weak coupling in 
the BGS regime. Here, in contrast, we investigate the 
intermediate to strongly coupled BCS-BEC crossover 
regime coupled with magnetic disorder, a problem that 
has not been explored so far, and which is now of very 
direct experimental relevance to atomic gases. We show 
that within Bogliubov-de Gennes mean field theory spin- 
dependent disorder exhibits qualitatively different fea- 
tures from the conventional symmetric case where both 
species feel the random potential. In particular, the gap 
and order parameter are driven to zero with increasing 
disorder, whereas they remain nonzero, and can even in- 
crease, when both species see the same energy landscape. 
This destruction of pairing has also been seen in BdG 
studies in the presence of spin-independent order and a 
uniform Zeeman field [9]. 

Spin dependent lattices can be achieved, for instance, 
by utilizing two atomic hyperfine transitions, Dl and 
D2 (ground to excited states) of the fermionic alkali 
atom ^^K, which can be individually excited by right- 
and left-circular polarized light respectively (according 
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to selection rules). By creating a lattice with inde- 
pendently tunable components of right- and left-circular 
polarized standing- wave light, a spin-dependent lattice 
is realized [l7|. The two fermionic "spin" species in 



this case are the two hyperfine states of the degenerate 
ground-state manifold split by a small uniform magnetic 
field. On the other hand, disorder can be superimposed 
on the periodic optical lattices created by interfering 
counter-propagating lasers by means fine-grained speckle 
fields 0, [HI. Accordingly, in principle, spin-dependent 
disorder can be realized by superimposing on the lattice 
separate laser beams that have been passed through a 
diffuser and focused at the lattice location. A suitable 
diffuser would be able to modulate the lattice intensity 
at each lattice site in a random way, with the randomness 
uncorrelated on the scale of lattice constants 18, loj . 
Model and Methodology: Our starting point is the at- 
tractive fermion Hubbard Hamiltonian with spin depen- 
dent disorder. 
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Here Cj^(cj^) are fermion creation (destruction) operators 
at site i for fermionic species (spin or hyperfine level) 
cr |. t is the hopping amplitude between near neigh- 
bor sites (ij) on a square lattice, t = 1 is the energy 
unit. —\U\ is an on-site attraction, and fia- is the chemi- 
cal potential, which in general will depend on the spin 
index a in order to maintain equal populations when 
the randomness is different for the two species. Cia is 
a local site energy which, in the conventional "Anderson- 
Hubbard" Hamiltonian, is chosen to be independent of 
fermion species a. Here we focus on the situation where 
= and is uniformly distributed on [— V, F]. 

In the clean limit, the square lattice attractive Hub- 
bard model is a supersolid (simultaneous long-range 
charge density wave and s-wave pair correlations) at 
T = and density p = I (half-filling). When doped 
to p 7^ 1, the degeneracy is broken, and the attrac- 
tive Hubbard Hamiltonian exhibits a finite temperature 
Berezinskii-Kosterlitz-Thouless transition to a supercon- 
ducting phase with Tc on the order of O.lt [20|. In the 
limit of strong attraction, this model crosses over to the 
Bose-Hubbard model of tightly bound pairs [21 1. 

Anderson's theorem provides a first insight into the 
effect of randomness on the clean attractive Hubbard 
model, suggesting that, at least for weak disorder, super- 
conductivity will survive. A number of numerical studies 
have quantitatively explored this model [22^. Here we 
will focus our comparison on the Bogoliubov-de Gennes 
(BdG) treatment of Ghosal et a/.0 since that is also the 
methodology employed here. 

The interaction term in H can be decoupled in different 



(charge, pairing, spin) channels. Since U < we focus 
on pairing and write, 

(ij),cr ia 



(2) 



where fiia- = Mcr + |^| (^i,-cr)- ^eff is diagonalized via the 
Bogoliubov transformation. 
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In the clean system the eigenfunctions Uin and Vin are 
plane waves. In the presence of disorder they are ob- 
tained by numerical diagonalization. The local order pa- 
rameter and density are determined self-consistently. 



Ai = -\U\{c,^c,^) = -\U\J2f{En)^ 



(4) 
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(where / is the Fermi function), as are the chemical po- 
tentials required to achieve the desired density |3]. These 
self-consistency conditions are equivalent to minimizing 
the free energy. The correctness and efficiency of our 
codes have been checked by using different initial config- 
urations of particle density and local pairing amplitude 
and checking for convergence. As expected, we found 
more iterations are needed for convergence near phase 
transitions. We define a spatially averaged order param- 
eter, Aop, from Ai. The BdG spectrum can be used 
to determine the energy gap, £^gap, which is the lowest 
eigenvalue above the chemical potential. The spectrum 
and eigenfunctions also detemine the density of states. 
Note that unlike the spin-independent case the eigenval- 
ues do not come in ± pairs and the distances of closest 
eigenvalues below and above the chemical potential are 
in general different. 

The key conclusions of the BdG treatment of Ghosal 
et ai 0] for the usual spin-symmetric case = are 
as follows: As disorder is increased the energy gap and 
lattice-averaged order parameter decrease, but never go 
to zero. In particular, even though the distribution of 
local pairing amplitudes has significant weight near zero, 
a finite spectral gap persists owing to remnant supercon- 
ducting islands and, importantly, significant overlap of 
the low energy excited states with these islands. Ander- 
son's theorem concerning the survival of superconductiv- 
ity is obtained under two assumptions: pairing of exact 
eigenstates and a further assumption that the kernel in 
the gap equation is spatially uniform. It is possible to 
generalize the calculation within pairing of exact eigen- 
states to allow for spatial structure in the pairing am- 
plitude. BdG calculations go beyond pairing of exact 



3 





/t (Tst-disorder) 

gap 

^A^p/t (Ti-disorder) 

/t (T-disorder) 
V gap " 
^^*.^^A^p/t (T-disorder)^ 






\ 

\ 

' — ' — 





0.5 



1.5 
V/t 



2.5 



FIG. 1. The lattice- averaged order parameter Aop and the en- 
ergy gap ^gap are shown as functions of the disorder strength 
V. In the spin- symmetric case considered in 0], these quan- 
tities never go to zero even at large V. In contrast, when the 
disorder is applied only to one spin species, sharp transitions 
are observed. The small "tail" in Eg^p goes to zero as the 
lattice size increases. N = 24 x 24, (n) = 0.875, U/t = -1.5. 



eigenstates and allow for a full treatment of the ampli- 
tude fluctuations in response to an underlying random 
potential. In order to see either the finite temperature 
phase transition from a superfluid to a non-superfluid 
state or the quantum phase transition from a superfluid 
to an insulator, thermal and quantum phase fluctuations 
must be included in the inhomogeneous BdG state. 

We will see that even at the level of pairing of exact 
eigenstates or BdG the situation is dramatically trans- 
formed by spin-dependent disorder. 
Results: Figure [T] shows the evolution of the energy gap 
£^gap and order parameter Aop with increasing random- 
ness V for density (n) = 0.875 and on-site attraction 
U = —1.5t. In the spin-symmetric case, confirming the 
results of , we find £^gap and Aop do not vanish. How- 
ever, in the spin-asymmetric case, when = 0, there 
are instead sharp transitions for both quantities. Inter- 
estingly, £^gap and Aop do not vanish simultaneously, but 
instead three phases are present: a superconductor char- 
acterized by Egg,p and Aop both nonzero at small disorder 
V, a gapless superconductor Igt [H, [1^ at intermediate 



y, and a E. 



gap 



phase at large V. 



The gapless superfluid has also been observed in other 
situations in which spin symmetry is broken, e.g. the 
mismatched Fermi surface considered in [25|. In such 
situations, paired regions of the Fermi sea coexist with 
regions occupied by only a single spin species. These un- 
paired pockets lead to a gapless superfluid phase between 
the normal and BCS regimes, much as occurs in Fig. [H 
Indeed, as we shall demonstrate later, in the attractive 
Anderson Hubbard Hamiltonian considered here, the real 
space distribution of Ai, indicates coexisting real space 
domains which are in direct analogy to the coexisting 
momentum space domains of [25| and related work. 




FIG. 2. The energy gap -E-gap, average order parameter Aop 
(left 2/-axis) and chemical potentials /i^ (right ?/-axis) are 
shown as functions of disorder strength V. Panels (a,b,c) 
show results for the case when —V<ei^<V and ai = 
for three values, U/t — —2, —3, —4. Panel (d) shows a case 
when both da are uniformly distributed on [— V^, V^] but are 
chosen independently, for U/t = —2. The qualitative behav- 
ior for these two types of spin-asymmetric disorder is similar: 
There are two transitions, first from gapped to gapless su- 
perfluid and then to a phase in which Aop = as well. As 
expected, /j.a is spin-independent in the situation of panel (d) 
when both species see disorder of the same overall strength 
(although distinct realizations). = 24 x 24, (n) = 0.875. 
The results are averaged over 10 realizations of randomness. 
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FIG. 3. Left panel: The distribution of the local order pa- 
rameter Ai. Right panel: The first excited state wave func- 
tion. Here U = —3t and = 4t. The average order parame- 
ter is non-zero, but because the first excited state has weight 
in the region where the local values are zero there is no energy 
gap. Hence the system is a gap less superfluid. 



Figures 2(abc) provide further details for the evolu- 
tion of £^gap and Aop for three interaction strengths, 
U/t = —2.0, —3.0 and —4.0. The transition points move 
to larger disorder strength, as expected, as the attractive 
interaction is increased. At all /7/t, there is a residual 
nonzero value of £^gap which scales to zero with increasing 
system size. In fact, the transition point from gapped to 
gapless superconductor is best signalled by a bifurcation 
of the chemical potentials fi^ 7^ fi^ required to maintain 
a spin-balanced population = p^. 

In order to determine whether such transitions are 
generic to the breaking of spin symmetry in the disor- 
der we consider, we show in Fig. 2d a case where both 
species see disorder drawn from the same distribution 
[— V] but with different realizations. That is, spin t 
and spin | fermions see different disordered landscapes. 
The basic features are the same as the case when only 
one species sees the disorder (Figs. 2abc): there are two 
separate transitions at which £^gap and Aop ^ 0. As 
expected, the chemical potentials do not bifurcate in the 
case of Fig. 2d. Nevertheless, there is still a signature 
of the vanishing of £^gap in the scatter of the chemical 
potential data. 

To understand the physics of the gapless superfluid 
we show, in Fig. 3, the spatial distribution of the lo- 
cal order parameter Ai (left) and the first excited state 
wave function (right). The system is highly inhomoge- 
neous and, although many sites have nonzero Ai, so that 
the average order parameter is substantial, there are also 
large domains where the local order parameter vanishes. 
The first excited state lives predominantly on the latter 
set of sites, leading to a vanishing energy gap. In the 
gapped superfluid phase at weaker disorder, ^* substan- 
tially overlaps regions of non-zero order parameter. Sim- 
ilarly, in the canonical case of spin-independent disorder 
the sites at which the first excited state has non-zero am- 
plitude coincide with those of non-zero order parameter. 

We can also infer the phase from the density of states 
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FIG. 4. Density of states at [/ = — 3t for (top to bottom) 
14 = 1.0t,2.0t and 3.0t. The gap closes for = 3.0t. See 
also Fig. 2(b). 

N{uj) obtained as a histogram of the BdG eigenvalues. 
Note again that the breaking of spin symmetry destroys 
the usual appearance of eigenvalues in ±£^n pairs. Figure 
4 shows the density of states 

n,i 

Nd^) = ^^\vin?5{co + E„) (5) 

n,i 

for U = —3t and a range of disorder strengths. The 
gap is nonzero for = IM but has closed by the time 

= 3.0t, in agreement with Fig. 2(b). 
Discussion: The effects of spin-dependent disorder on the 
superconducting phase in the attractive Hubbard Hamil- 
tonian differ fundamentally from conventional chemical 
potential disorder. Sharp transitions at which the energy 
gap and average order parameter vanish can clearly be 
seen as the randomness increases. This can be attributed 
at the fundamental level to the breaking of time reversal 
symmetry, which can arise in the solid state by the spin- 
dependent scattering produced by magnetic impurities. 
In fact, the spin-dependent site energy ei^ni^ + ^ii^ii in 
Eq.[2l can be thought of as a combination of a local chem- 
ical potential ^(ci^ + ei^) {n[^ + ni^) and a local Zeeman 
field, ^(cit — Cii) {rii^ — riii). The Zeeman term breaks 
time reversal symmetry, providing one way to interpret 
the results presented here. The occurrence of sharp tran- 
sitions can also be understood in analogy with exotic 
paired states in other situations like Fulde-Ferrell-Larkin- 
Ovchinnikov(26l [27! with mismatched Fermi surfaces (iij. 
In these systems, electrons of one species with a particu- 
lar momentum cannot find partners, and hence a super- 
fluid and unpaired electrons coexist. An analogous phe- 
nomenon occurs in real space in the Hamiltonian consid- 
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FIG. 5. The behavior of the energy gap and the order pa- 
rameter for the same disorder strength for the case when the 
local order parameter is forced to be uniform to conform to 
the conditions for which the classic Abrikosov-Gorkov the- 
ory [3l was developed and its comparison with a spatially 
non-uniform pairing amplitude. Notice that a finite order pa- 
rameter can exist up to much higher disorder when the pairing 
amplitude is allowed to be non- uniform. N = 24 x 24, U /t — 
-2, (n) = 0.875. 



ered in this paper: the spin-dependent randomness pro- 
vided regions of the lattice where fermions cannot find 
partners, while in other regions pairing can continue. 
Comparison with Abrikosov-Gorkov theory: It is useful 
to compare our results in Figs. 2 and 3 against the early 
work by Abrikosov and Gorkov (AG) [16] where they first 
predicted the possibility of gapless superconductivity for 
spin-dependent scattering, a phase with finite order pa- 
rameter but no gap. Their result was obtained within 
perturbation theory in a regime of weak disorder assum- 
ing that the local pairing amplitude is homogeneous. 

In Fig. 5 we compare our results for the energy gap and 
order parameter obtained from the inhomogeneous BdG 
equations as discussed above against a restricted calcu- 
lation where we enforce a spatially uniform local pairing 
amplitude A = {1/N) A(ri) by finding a self consis- 
tent solution for the average pairing amplitude A. The 
significant aspects of these results are: (1) When the local 
pairing amplitude is forced to be uniform, there is a small 
region in disorder strength where a gapless phase exists, 
in agreement with AG theory 



16|. (2) When the local 



pairing amplitude is allowed to vary spatially, a much 
larger region of gapless superconductivity opens up. The 
gap is not affected much by the local pairing amplitude 
variations but the disorder strength up to which a finite 
order parameter exists as deduced from long range su- 
perconducting correlations is greatly enhanced. 

In that respect the gapless superconducting phase that 
is reported here is different in nature from that in the 
earlier work because inhomogeneity of the local order pa- 
rameter plays a vital role in enhancing the order param- 
eter. 

Conclusion: A key goal of optical lattice emulation is to 



observe Quantum Phase Transitions (QPTs) in ultra-cold 
atomic gases, analogous to those believed to underlie con- 
densed matter phenomena such as high temperature su- 
perconductivity. Our work suggests that, for exploring 
QPTs associated with the interplay of attractive inter- 
actions and disorder, the sorts of spin- dependent lattices 



25| might be especially 



achieved experimentally 
useful. 

Several fundamental conceptual questions remain to be 
addressed by these experiments, and by more precise the- 
oretical approaches such as Quantum Monte Carlo. One 
is the construction of a unified phase diagram in which 
spin-dependence arising from the scattering, spin depen- 
dent disorder considered here, and from imbalanced pop- 
ulations, FFLO physics, are both included. Specific is- 
sues include the conditions for the existence of gapless su- 
perfluidity and whether an inhomogeneous order param- 
eter is essential for its formation. Another outstanding 
question concerns transport properties: Starting from the 
weakly interacting limit, in which one species feels dis- 
order and the other does not, the system is 'half-metal' 
with distinct localized and itinerant components. As an 
attractive interaction U is turned on, the delocalized par- 
ticles will feel an induced randomness as they interact 
with the localized species. Is there a critical U separat- 
ing distinct half-metallic phases from phases in which the 
two components have identical transport characteristics? 

We acknowledge financial support from ARO Award 
W911NF0710576 with funds from the DARPA OLE 
Program, the CNRS-UC Davis EPOCAL LIA joint re- 
search grant, the Department of Energy under DE-FG52- 
09NA29464, and useful input from B. Brummels. 



[1] 
[2] 

[3: 

[4] 
[5] 
[6] 

[7] 

[8] 
[9] 

[lo: 
[11 

[12 

[13 

[14 
[15 



A. Cho, Science 320, 313 (2008). 

P.W. Anderson, Phys.Rev. 109, 1492 (1958). 

P.W. Anderson, D. C. Licciardello, and T. V. Ramakr- 

ishnan, Phys. Rev. Lett. 42, 673 (1979). 

A. Aspect and M. Inguscio, Physics Today 62, 30 (2009). 

M. White et a/., Phys. Rev. Lett. 102, 055301 (2009). 

L. Fallani et al, Phys. Rev. Lett. 98, 130404 (2007). 

A. Ghosal, M. Randeria, and N. Trivedi, Phys. Rev. Lett. 

81, 3940 (1998); Phys. Rev. B65, 014501 (2001). 

K. Bouadim, Y.L. Loh, M. Randeria Nature Physics 7, 

884 (2011). 

Y. Dubi, Y. Meir and Y. Avishai, Phys. Rev. B78, 
024502 (2008). 

M.P.A. Fisher, PB. Weichman, G. Grinstein and D.S. 
Fisher, Phys. Rev. B40, 546 (1989). 
W. Krauth, N. Trivedi and D. Geperley, Phys. Rev. Lett. 
67, 2307 (1991). 

R.T. Scalettar, G.G. Batrouni, and G.T. Zimanyi, 
Phys. Rev. Lett. 66, 3144 (1991). 

S.G. Soyler, M. Kiselev, N. V. Prokofev, and B. V. Svis- 
tunov, Phys. Rev. Lett. 107, 185301 (2011) 
S. Iyer, D. Pekker, G. Refael, arXiv: 11 10.3338] 
M. Randeria, Nature Physics 6, 561 (2010). 



6 



[16] A.A. Abrikosov and L.P. Gor'kov, Soviet Physics JETP 

12 1243 (1961). 
[17] W.V. Liu, F. Wilczek, and P. Zoller, Phys. Rev. A70, 

033603 (2004). 

[18] O. Mandel et al, Nature 425, 937 (2003); Phys. Rev. 

Lett. 91, 010407 (2003). 
[19] D. McKay and B. DeMarco, New Journal of Physics 12, 

055013 (2010). 

[20] S. Robaszkiewicz, R. Micnas, and K.A. Chao, 
Phys. Rev. B24, 1579 (1981); R.T. Scalettar et al, 
Phys. Rev. Lett. 62, 1407 (1989). 

[21] M.P.A. Fisher et al, Phys. Rev. B40, 546 (1989). 



[22] N. Trivedi, R.T. Scalettar, and M. Randeria, 

Phys. Rev. B54, 3756 (1996). 
[23] P.G. de Gennes, Superconductivity in metals and alloys 

(Benjamin, New York, 1966). 
[24] A. Lamacraft and B.D. Simons, Phys. Rev. Lett. 85, 4783 

(2000). 

[25] A.E. Feiguin and M.P.A. Fisher, Phys. Rev. Lett. 103, 
025303 (2009). 

[26] P. Fulde and A. Ferreh, Phys. Rev. 135, A550 (1964). 
[27] A. Larkin and Y.N. Ovchinnokov, 

Zh. Eksp. Teor. Fiz. 47, 1136 (1964) [Sov. Phys. JETP 

20, 762 (1965)]. 



